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Abstract 

We analyse second moments i?2 of like-sign pion pairs in the two-dimensional {cil,(1t) and 
three-dimensional (qoiQSiQL) decompositions of the three-momentum difference. Conven- 
tional fit parametrisations such as gaussian, exponential, power-law and Edgeworth fail mis- 
erably, while more elaborate ones such as Levy do well but fail to yield a unique best-fit 
solution. A two-component model using a hard cut to separate small- and large-scale parts 
appears possible but not compelling. In all cases, the data exhibits a strong and hitherto 
unexplained peak at small momentum differences which exceeds current fits. 

Classification: 13.85.Hd, 13.87.Fh, 13.85.-t, 25.75.Gz 
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The UAl experiment, having completed data-taking at the CERN SPS in the late eighties, 
continues to be relevant and interesting. The current concentrated effort at RHIC to quan- 
tify and understand ultrarelativistic nuclear collisions relies extensively on comparisons with 
baseline scenarios constructed from the corresponding "trivial" hadron-hadron sample. Current 
experimental energies of 200 AGeV at RHIC are still below those available to UAl by a factor 
three, so that UAl results may also provide a window on possible energy dependencies of current 
investigations. 

In this contribution, we provide preliminary results on HBT analysis mainly in terms of the 
two-dimensional decomposition of the three-momentum difference, defining in the usual way 
Ql = IqlI = Kcf-z) z\, with z the beam direction, and qt = |q — QlI- Brief reference is also made 
to the three-dimensional Bertsch-Pratt case, with similar results and issues arising. 



Like-sign (LS) pion pairs from approximately 2.45 million minimum-bias events measured by 
the UAl central detector were analysed. This represents a twofold increase in statistics over 
Ref. P and a 15-fold increase compared to earlier UAl HBT analyses 2t^^. Standard cuts 
^ were applied, including single-track cuts p± > 0.15 GeV/c, \y\ < 3 and, to avoid acceptance 
problems, 45° < |0| < 135°. The sample contains an estimated 15% contamination of charged 
kaons. 



The most important among the standard pair cuts is the "ghost cut" which eliminates spurious 
"split track" LS pairs within a narrow cone but many real LS pairs also. A correction fac- 
tor compensating for this was determined by passing unlike-sign (US) pairs through the same 

^Contribution to the Workshop on Particle Correlations and Femtoscopy, Kromefiz, Czech Republic, 15-17 
August 2005. To appear in: Proc. 35th International Symposium on Multiparticle Dynamics, Kromei'iz, 9-15 
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ghost-cut algorithm. Such correction factors were determined for each (q^, qx) bin^ and charged- 
multiphcity subsample; for some {qL,QT) bins, this correction factor ranges up to 1.7 or even 1.9 
for low-multiphcity subsamples. 

We corrected for Coulomb repulsion by parametrising the Bowler Coulomb correction in the 
invariant momentum difference Q = \/—{pi — V'l)'^ with an exponentially damped Gamov 
factor G(Q) 6., with a best-fit value Qeff = 0.173 ± 0.001 GeV/c, 



i^coui(Q) = 1 + \G{Q) - 1] exp(-Q/Qeff) 



(1) 



The reference sample was formed by randomly combining LS tracks taken from pools of events 
in the same subsample of event charged multiplicity as the sibling event currently being 
analysed. Note that the reference for fixed-multiplicity subsamples is not the poisson distribution 
but the multinomial, whose second moment is pf'"l'(q|iV) = (1 - iV-i)pi8)pi(q| A^), so that 
the appropriate normalised moment is ^ [7| 



i?2(q) 



(2) 



pi! 
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The (ql,<7t) bins for ghost corrections use qh measured in the detector rest system (pp CMS), while all HBT 
quantities are measured in the LCMS. 
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FIGURE 1. Upper panels: R2{qL,qT) data and best fits, shown for the first shces (left to 
right) qr = 0.00-0.02, 0.02-0.04, 0.04-0.06 and 0.06-0.08 GeV/c. Lower panels: R2{qL,qT) 
and the same fits for corresponding fixed-q^ slices. Solid lines: Gauss/Edgeworth fits; dashed: 
exponential; dotted: exponential with cross term; dash-dotted: Levy. The bin {qL,qT) < 
(0.02, 0.02) GeV/c is omitted from all fits and plots. 



The ghost- and coulomb-corrected normalised moment R2{qL,QT) is shown in Fig. 1, together 
with fits to parametrisations R2 = 7[1 + A|S'i2p], with |>S'i2|^ parametrised respectively as 
exp(— — R^q^ — 2RLTqLqT) (gauss with cross term), exp{—RLqL — Rrqr) (exponential), 



eyiY>{— Riqi— RrqT—'^RLT ^/qUlT) (exponential with cross term) and i?2 = 7 
(product power law). Note that all results shown are preliminary. 



1 + (i^L^L)""^ {RxqT) 



-ax 



It is immediately apparent that none of these fits reproduces the data, with x^/NDF ranging 
from 3 to 9. A parametrisation based on an Edgeworth expansion [11111) 



\S- 



12 



= exp{-Rlql-R^q^) [] [l + K^^a H^{V2 Rdqd)/2i 

d=L,T 



(3) 



(with K4^rf the fourth-order cumulant in q^ and the corresponding hermite polynomial) fares 
no better as best-fit values for ^4^^ in both directions turn out to be negligible. We omit the 
third-order terms in k^^l and k^^t as they are antisymmetric in qa- 
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FIGURE 2: Comparing the shapes of R2{qL,qT) at intermediate scales. Panels show (a)) 
UAl data, (b) gauss, (c) exponential, (d) power-law, (e) Levy, (f) exponential with cross 
term fits. The gauss-type fits have the right shape, but end up far below the data peak 
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at small {qL,qT)- While the simple exponential and product power-law may approximate 
the peak reasonably, both fail miserably when shapes at intermediate scales are considered. 
Levy does well but no unique best fit can be found. Shape-wise, the exponential with cross 
term (f) appears to come out on top. All plots are truncated at R2 < 1.9 in order to bring 
out structure at intermediate scales. 



A Levy-based parametrisation jlOj . 

\Su\' = eM-Rhl-RWTr^^ (4) 

yields better results in reproducing the strong peak observed in the data; however, four of the 
five fit parameters, viz. A, Rl, Rt and a, are strongly correlated so that no unique best fit 
can be achieved. One example of many equivalent "best" fits is shown in Fig. 1. Omitting a 
second small- ( g^,, qt) bin from the Levy fit renders the fit even more unstable (the innermost 
bin {qltIt) < (0.02,0.02) GeV/c is excluded from all analysis as a matter of course). This is 
hardly surprising, as the four abovementioned parameters collectively depend strongly on the 
exact shape of the peak in the very small (qljQt) region — the very region that experimental 
measurement struggles to resolve. 

Interesting, nonetheless, is the observation that fits fail for different reasons: As shown in Fig. 2, 
the gauss and Edgeworth fits reproduce the shape of R2 at intermediate (qljQt) rather well 
but lacks the strong peak exhibited in the data,^ while the exponential and product power-law 
parametrisations are more peaked but fail to describe the shape at intermediate {qljQt)- Judg- 
ing by shape alone, the exponential with a 2Rlt^/(1l Qt cross term does best (Fig. 2(f)). 

Either way, it is clear that the superiority of the power-law fit to the one-dimensional invari- 
ant four-momentum moment R2{Q) over Gauss and exponential fits |^ is not repeated in the 
R2{qL,qT) case due to the hyperbolic shape of our "product power-law" parametrisation^ com- 
pared to the elliptic form of the R2{ql,Qt) data shown in Fig. 2(a). The analysis of shapes 
at intermediate scales appears to provide valuable information. Shape analysis in the form of 
Refs. may help to quantify these qualitative observations. 

The failure of conventional parametrisations to reproduce the elusive peak at small (qljQt) 
suggests that there may be two scales in the system. Borrowing the (strictly speaking inappro- 
priate) terms "core" and "halo" from the literature we try a two-scale "butcher's model", 
similar to simpler precursors in Refs. |141 115j . 

R2{qL, qr) = 1 l + Ac e^p{-Rlcql - R^cIt) + e^pi-Rl^ql - R^h(1t) > (5) 

and proceed as follows: First, we fit only bins with momentum differences larger than a hard 
cutoff, ((?L,9r) > (^cut, ^cut), to the core gaussian. The resulting best-fit is subtracted from all 
data. The remaining halo "data" with {qL,qT) < ('Zcuti Qcut) is then fit with the halo gaussian, 
and the resulting "core" and "halo" best fits combined. Data and fit histograms at various steps 
of this procedure are shown in Fig. 3. 



''Edgeworth fits are indistinguishable from the normal gaussian fits throughout this analysis. 

■^Other forms, such as a "sum power law" _R2 = 7 [l -I- {Rl <?!-) + {Rt ^/t) remain to be tested. 
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FIGURE 3: Two-scale butcher's model. Of the i?2(<ZLi'?T) data (left upper panel) only 
bins with (qLiQT) > (0.16,0.16) GeV/c are fit to the "core" gaussian of Eq. ^ as shown 
in the right upper panel. The result is subtracted from the data to reveal the "halo" (right 
lower) which is then fit to the halo part of Eq. (0). The combined core-halo fit is shown in 
the left lower panel. 
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FIGURE 4: Dependence on gcut- (a) Total for both core and halo fits as a function of 
Qcut- The two lowest values correspond to <7cut = 0.16 and 0.18 GeV/c, giving x^/NDF = 
2.28. Panels (b)-(e): Dependence of best-fit parameter values on gcut- Squares and circles 
represent core and halo fit parameters respectively. A clear separation of scales is observed 
as assumed in the model. The halo continues to have a chaoticity parameter Xh that exceeds 
its theoretical limit of 1. 



In Fig. 4, the dependence of the two-scale model on gcut is tested. As shown in Fig. 4(a), the 
joint for both fits as a function of ^cut is found to have a well-defined minimum for qcut = 
160-180 MeV/c. The best estimates for R^, Rt, A and Rl/Rt are the two filled points in 
Fig. 4(b)-(e) corresponding to the two best q^nt values. Averaging these two numbers, we esti- 
mate Rlc = 0.82 ± 0.02 fm, Rtc = 0.49 ± 0.02 fm, Rlh = 2.45 ± 0.03 fm, Rth = 2.29 ± 0.05 
fm, Ac = 0.26 lb 0.01 and Xh = 1-17 it 0.02, signalling a prolate "core" and a roughly spherical 
"halo" . The clear separation between the sizes of the core and halo radii a posteriori support the 
assumption of the presence of two scales, i.e. signal that the two-scale model is consistent. We 
note that Xh continues to exceed the theoretical limit of 1.00, albeit not as strongly as the huge 
intercept i?2(0,0) > 2.7 seen in the data itself. The joint best x^/NDF = 2.28 for the two-scale 
model is still rather large, however, so that all numbers and conclusions should he treated with 
caution. 

Turning briefly to the more common three-dimensional Bertsch-Pratt representation ) we find 
that the data once again has a strong peak at small {qo,Qs,QL), and that the simple gaussian 
parametrisation fails completely. In Fig. 5, we show the result of fitting the Levy and Edgeworth 
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parametrisations, 



\s,,f = eM-RUh-RWs-Rhlr^', (6) 

= eM-RoQo-RsQs-Rhl) n [l + ^^4,dH4(.V2Rdqd)/24\, (7) 

d=0,S,L 

to R2{QOyQSTQL), with slices plotted along the three axes of the three-dimensional space. Even 
these parametrisations appear not to describe the data well for small |q|, and in this threedi- 
mensional case the discrepancies are spread more widely than in the {qL,QT) case. Also, while 
the Levy fit does appear to do reproduce the data to a reasonable degree, it again suffers from 
strong inherent correlations between the parameters A, Ro,Rsj Rl and a. Such correlation again 
implies that no unique minimum for the exists and thereby no unique best-fit values for its 
parameters either. 
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FIGURE 5: Second moment R2{qo , Qs i Ql) of three-dimensional momentum difference 
decomposition, plotted along the three axes qd, with the other two variables {qe,qf) < 
(0.02,0.02) GeV/c. Solid Hne: Levy fit; dashed Hne: Edgeworth. 



It is, of course, always possible to try even more elaborate parametrisations, for example using 
separate a-exponents for each of the {qo,Qs,QL) directions. While a better fit might theoretically 
be achieved, this will invariably come at the price of even more highly correlated parameters. 
The data cannot distinguish between many combinations of "best" values for such parameters. 

We note that, even within the arguably artificial method of the two-scale model, the strong peak 
in the data is not reproduced, resulting in the quoted low confidence level. Nevertheless, it may 
provide a useful hint that "something else is going on" , be it the influence of jets, clustering 
effects or some other unknown factor. 



We stress that is is unlikely that the strong peak seen in our data at small momentum differ- 
ences, over and above simple gaussian or exponential parametrisations, is due to bias. First, the 
peak persists even without the Coulomb or split-track corrections. Second, we believe that the 
same peak may well be responsible for the power-law form and higher-order effects seen in our 
earlier work with one-dimensional distributions Third, the presence of unidentified kaons 
and protons in the sample imply that the real peak should exceed the one shown here, so that 
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the present results may be lower limits. It should also be noted that a number of other hadronic 
experiments have previously seen significant deviation from gaussian behaviour, especially at 
small I q| jl71 118j . The challenge is clearly now to find a convincing physical cause and explana- 
tion. 
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